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Abstract 

We perform a coupled-channel calculation of the H dibaryon within a chiral constituent quark 
model. The problem is solved within a quark model constrained by the experimental data of 
strangeness -1 and -2 two-baryon systems. We examine in detail the role played by the different 
contributions of the interacting potential as well as the number of coupled channels considered. 
Special attention has been payed to the parameter dependence, flavor symmetry breaking and spa- 
tial configurations. The value extracted for the binding energy of the H dibaryon, being compatible 
with the restrictions imposed by the Nagara event, falls within a plausible extrapolation of recent 
lattice QCD results. 

PACS numbers: 12.39.Jh,12.39.Pn,14.20.Jn,14.20.Pt 
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I. INTRODUCTION 



Very recently the H dibaryon was put back on the agenda by lattice QCD calculations 
of different collaborations, NPLQCD [lj and HAL QCD j2j, providing evidence for a bound 
state at pion masses larger than the physical ones. It was first proposed thirty years ago ^ 
as a spin and flavor singlet composed of six quarks (uuddss). Such a proposal emerged on 
the context of the MIT bag model [4j . When applied to all six-quark systems it predicted 
the existence of only one stable dihyperon with J p = + and a mass of 2150 MeV, 81 MeV 
below the AA threshold, therefore strongly bound. Moreover, being the lightest particle of a 
two-baryon system with strangeness (S) -2, it would be stable against the strong interaction 
and would necessarily decay weakly. It was assumed to be a single hadron made of six quarks 
squeezed in a small region and not a two-baryon state bound in an S wave like the deuteron. 
If this were the case it would open the door to the exotic states, i.e., hadrons that do not fit 
in the standard qq or qqq configurations. Otherwise, the fact that the H dibaryon is the most 
promising candidate to be the second dibaryonic state after the deuteron has generated a 
lot of expectation. 

On the experimental side, there are very few data in the S = —2 sector, coming from the 
inelastic S~p — > AA cross section at a lab momentum of around 500 MeV/c, and from the 
elastic H~p — > S~p and inelastic — y S°n cross sections for lab momenta in the range of 
500 - 600 MeV/c - Thus, the relevant information we have is indirect and comes from 
double-A hypernuclei. Their binding energies, -Baa ; provide upper limits for that of the H 
dibaryon, i.e., B_u < -Baa- The first hypernuclear events are quite old and admit several 
interpretations In 2001 it was reported the so-called Nagara event [11], interpreted 



uniquely as the sequential decay of ^He emitted from a H~-hyperon nuclear capture at rest. 
The mass and the values of -Baa and of the AA interaction energy, AB AA , were determined 
without ambiguities. The small value of AB^a suggested an attraction weaker than the 
one previously estimated. It also gave the most stringent constraint to the mass of the 
H dibaryon to date, i.e., Mh > 2223.7 MeV at a 90% confidence level. It took almost 
one decade, but four more double-A hypernuclear events were reported, from KEK E176 



and E373 experiments [12] , still with preliminary results. All the details are summarized in 



Table [U The future E07 experiment from J-PARC [12| is expected to improve our knowledge 



TABLE I: Reported double hypernuclear events. 



Event 


Nuclide B AA (MeV) AB AA (MeV) 


1963 


i°ABe 


17.7 


± 


0.4 


4.3 


± 


0.4 


1966 


AAHe 


10.9 


± 


0.5 


4.7 


± 


1.0 


1991 


13 R 
AA a 


27.5 


± 


0.7 


4.8 


± 


0.7 


NAGARA 


AA^ 


7.13 


± 


0.87 


1.0 


± 


0.2 


MIRAGE 


AAHe 


10.06 


± 


1.72 


3.82 


± 


1.72 


DEMACHIYANAGI 


i° A Be 


11.90 


± 


0.13 


-1.52 


± 


0.15 


HIDA 


AABe 


20.49 


± 


1.15 


2.27 


± 


1.23 




12 Re 

AA 


22.23 


± 


1.15 








E176 


13 T3 p 

AA 


23.3 


± 


0.7 


0.6 


± 


0.8 
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on the S =-2 sector, giving ten times more events. 

From the theoretical point of view, many approaches have been performed, their predic- 
tions for the binding of the H spreading over a wide range of energies [l3|. Recent lattice 
results produced by the NPLQCD Q and HAL QCD Q Collaborations found a bound H 
dibaryon for non-physical values of the pion mass (m n = 837 MeV and = 670 — > 1010 
MeV, respectively). When performing quadratic and linear extrapolations to the physical 
point [llj], a bound dibaryon (around 7 MeV) and a H at threshold, respectively, are pre- 
dicted. Also presented in Ref. [lij are preliminary results for m n = 230 MeV, much closer 
to the physical pion mass, pointing to a H dibaryon at threshold, as also experimentally 
suggested by the enhancement of the AA production near threshold found in Ref. Ts| . 



In this work we present the first approach to the H dibaryon within a constitue nt q uark 



model constrained by the experimental data of the S — -1 and -2 cross sections [16j, [17 
To make our results more robust and significative, we scrutinize all channels in the S = -2 
sector, paying due attention to the H dibaryon, (T,S) = (0,0). The paper is organized as 
follows. In Section II we provide a brief description of the constituent quark model and the 
formalism to study the coupled-channel problem. In Section III we present and analyze our 
results. Finally, we summarize our conclusions in Section IV. 



II. THE CHIRAL CONSTITUENT QUARK MODEL 

The baryon-baryon interactions needed for the study of the H dibaryon are computed 
from a chiral constituent quark model (CCQM) [18]. Baryons are described as clusters of 
three effective constituent quarks, their mass coming from the spontaneous chiral symmetry 
breaking. The CCQM was first applied to the study of the nonstrange SU(2) xSU(2) sector, 
describing the baryon spectroscopy and the NN interaction in a consistent manner (l8| . 
Such a success was due to the choice of the adequate mechanisms in the description of 
the quark-quark-meson interaction: quark antisymmetry plus a perturbative short-range 
force together with a microscopic non-perturbative chiral interaction at medium and long 
distances. They allow a unified treatment of the one-, two- and three-body systems, with 
a reduced and unique set of parameters. Later on, a generalization to SU(3)xSU(3) was 
done to perform a systematic and detailed analysis of the qq spectrum [l9[ . The interaction 
between quarks was given by (Model I from now on): 

V qq {f) = VconW + Vb GE (r) + V a {r) + V PSE {f) , (1) 

where Vcon is a confinement term that represents the nonperturbative aspects of QCD, Voce 
is the one-gluon exchange (OGE) potential, obtained through the nonrelativistic reduction 
of the quark-quark-gluon interaction diagram in QCD, V a is the one-sigma-exchange po- 
tential and Vpse stands for the chiral potential, associated to the exchange of pseudoscalar 
Goldstone bosons. It comprises one-pion, one-eta and one-kaon exchanges. When studying 
the S — —1 Ap cross section, it was noticed that a better approximation to the scalar inter- 
action was needed, considering the nonet of scalar mesons, that comprises a singlet (denoted 
by <7 ) and an octet that will be denoted by Vsce, giving thus rise to Model II: 

V w (f) = V CON (r) + V OGE (r) + V ao (r) + V SCE (f) + V PSE (f) , (2) 

where the relation between a and do is given by a = cos9 s uq + sin# s a 8 . The present model 



has been used in Ref. [16j to study two- and three-baryon systems with strangeness -1 
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giving a nice description of the hypertriton. It has also been used to study the strangeness 
-2 two-body scattering cross sections [17| . A reasonable fit to the experimental data of 
the elastic S~p and the inelastic S~p — > rPn and — > AA cross sections reported in 
Refs. 

was obtained. 

In order to derive the local BiB 2 — > B^B^ potentials from the basic qq interaction we use 
a Born-Oppenheimer approximation. Explicitly, the potential is calculated as follows, 

V Bl B 2 (LST)^B 3 B 4 (L'S'T)(R) = ClST T ( R ) ~ ^LS^i 00 ) > ( 3 ) 

where 



\*%k T &) i Ef<y=i w^o i mm, 

£lst ( r ) = i i ■ ( 4 ) 

In the last expression the quark coordinates are integrated out keeping R fixed, the resulting 
interaction being a function of the Bi — Bj relative distance. The wave function ^^b{R) 
for the two-baryon system is discussed in detail in Ref. (T^J. This formalism allows us to 
isolate different contributions and/or diagrams, making it easy to analyze the results. 

Once we have the two-baryon interactions, we switch to solve the two-body coupled- 
channel problem. Let us start from a physical system made of two baryons, B\ and B 2 
(Bi = A,N, S, E), with isospin, spin and parity quantum numbers (I)J P in a relative S 
state. They interact through a potential V that contains a tensor force. Then, in general, 
there is a coupling to the B\B 2 D— wave and to any other two baryon system (AA, NE, 
EE) that can couple to the same quantum numbers (I)J P . Thus, if we denote AA = D\, 
NE = D 2 and EE = D3, the Lippmann-Schwinger equation for the S = —2 BiB 2 scattering 
becomes 

POO 
_n o^O 



X 



,0,2' 

1=D k 
(fc=l,2,3) 

G^EipJtfyfWfaME) ,a,/3 = D 1 ,D 2 ,D 3 , (5) 

where t is the two-body scattering amplitude, j, i, and E are the angular momentum, isospin 
and energy of the system, respectively; £ a s a , £ 7 s 7 , and l$sp are the initial, intermediate, and 
final orbital angular momentum and spin, respectively, and p 7 is the relative momentum of 
the two-body system 7. The propagators G 7 (i?;p 7 ) are given by 

with 

E=^, (7) 
2/x 7 

where /i 7 is the reduced mass of the two-body system 7. For bound-state problems E < 
so that the singularity of the propagator is never touched and we can forget the ie in the 
denominator. If we make the change of variables 
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where b is a scale parameter, and the same for p a and pp, we can write Eq. (JSJ) as 



2 ^ 7 

l=D k 7 
(fc= 1,2,3) 

x (x a , x 7 ) G 7 (£; p 7 ) (* 7 , x,; £) . (9) 

We solve this equation by replacing the integral from —1 to 1 by a Gauss-Legendre quadra- 
ture which results in the set of linear equations 

N 

E E E<^ s w5^ (io) 

£~=0,2 m=l 
7=£>fc ' 

(fe=l,2,3) 

with 

2 



V 1 - X m / (1 - X r 



x ^ ; J V ^ n ,x m )G 7 (£;j, 7 J, (11) 

and where w m and the weights and abscissas of the Gauss-Legendre quadrature, 

while p lm is obtained by putting Eq. (|S]). If a bound state exists at an energy 

Eb, the determinant of the matrix M^J"' ~ <Sl (Eb) vanishes, i.e., \M aTj ji(EB)\ = 0. We 
took the scale parameter b of Eq. (jSD as 6 = 3 fm _1 and used a Gauss-Legendre quadrature 
with TV = 20 points. 



III. RESULTS 

The flavor singlet (T, S) = (0, 0) state in the S = — 2 sector comprises three coupled 
two-baryon channels, presenting rich coupling effects. If SU(3) were exact, one could write 

W = \/^ AA > + vf m) - vf • (12) 

The presence of three transition potentials, namely AA — NS, AA — EE and NE — EE, 
makes it necessary, to obtain the bound states of the system, to perform a coupled-channel 



TABLE II: Single channel binding energy (in MeV). 





B AA 


Bne 


-Bee 


Model I 


1.9 


32.7 


24.2 


Model II 




0.1 


0.1 



5 



[I 



'R(fln) 



o 



1 2 

R(fin) 



FIG. 1: Interacting potentials contributing to the H dibaryon in (a) Model I and (b) Model II. 



calculation where six different interactions contribute: three diagonal potentials (AA, A"S 
and EE) plus the three transitions just mentioned. 

Regarding the baryon-baryon interactions, let us only point out here a couple of remark- 
able features of interest for the H dibaryon study. First of all, the three diagonal potentials 
are attractive. The smaller attraction corresponds to the AA potential, whereas NS and 
EE show broader and deeper attractions, respectively. In order to know about the magni- 
tude of these diagonal interactions, the Fredholm determinant -that tells the binding energy 
of a system- can help. The binding energies of these three channels with respect to their 
own thresholds (still without any coupling), computed in Models I and II, are shown in 
Table [Til The six potentials contributing to the H dibaryon in Models I and II are dis- 
played in Figs. Ufa) and (b). From them we shall proceed to solve the coupled-channel 
problem by enlarging the Hilbert space progressively, from {AA} to {AA, A r S} and finally 
to {AA, V5,EE}. The thresholds of these three channels are respectively 2231 MeV, 2257 
MeV and 2381 MeV. 

In this way we are able to separate and to identify the modifications introduced by every 
baryon-baryon channel. The results, computed in Models I and II, are shown in Table HTT1 
When approximating the interaction by the AA potential, one finds the H dibaryon slightly 
above threshold in Model II. As expected, due to the attractive character of the remaining 
channels, the coupling increases the binding energy of the dihyperon by several MeV. The 
low sensitivity of this binding when going from Model I to Model II contrasts with the 

TABLE III: S = —2, (T, S) = (0, 0) binding energy (in MeV) in the one-, two- and three-channel 
approximations . 





B{AA} 


B{AA,NE.} 


-E>{AA,7VH,EE} 


Model I 


1.9 


8.7 


10.0 


Model II 




1.6 


7.0 




FIG. 2: OGE interaction in the diagonal channels contributing to the H. 



important modification of the Ap cross section [l6j]. Hence a detailed analysis of the different 
contributions to the binding may be helpful. For that purpose we shall restrict ourselves to 
Model II, as it is the one that properly describes the S = — 1 and —2 experimental cross 
sections [171 ] . 



A. Pieces of the interaction 

A more precise analysis can be done by studying the different components of the interaction 
separately. Let us start with the chromomagnetic interaction consequence of the OGE, the 
basic mechanism originally proposed to bind the H [Jjj] . In the work by Oka et al. [20] , using a 
one-gluon exchange plus a confining potential, a resonance 152 MeV above the AA threshold 
was obtained. The large difference with respect to the original work [jjj], that gave Bh = 81 
MeV, was a consequence from the fact that Jaffe [3} used the limit of flavor SU(3) symmetry 
and took for (uuddss) the short-range correlation coefficients as for ordinary hadrons. Oka 



et al. |20| showed that in the dilute (uuddss) system the strength of chromomagnetic effects 
is reduced as compared to ordinary hadrons. 

Later on, several works employing one-gluon exchange potentials have approached the 
study of the H dibaryon offering again a wide range of results. Ref. j2l[ took into account two 
possibilities: six-quark and two-cluster configurations. In the first one, even in the simplest 
case of OGE plus confinement, the H dibaryon was more than 300 MeV above threshold. In 
the second one, although the H was not yet bound, the energies of the resonance were much 
smaller. This was due to the repulsive character of the chromomagnetic interaction in the 
AA potential, not compensated by its attraction in the remaining channels. The very same 
feature was found by Nakamoto et al. [22|, where, moreover, the repulsion in AA was much 
stronger. On the other side, Stancu et al. j23|, by performing a study of the short range 
part of the Goldstone boson exchange interaction, lacking thus the OGE potential, found 
a resonance more than 800 MeV above threshold. Finally, in a chiral quark model study 



with six-quark (Os) 6 as well as two-baryon configurations done by Shimizu et al. [24|, three 
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TABLE IV: Character of the strangeness -2 two-baryon interaction in the one-, two- and three- 
channel approximations, for different quark-quark interactions. R indicates repulsion, WA weak 
attraction, A(N) indicates attraction, being N the binding energy in MeV. PSE stands for the 
pseudoscalar exchange, OGE for gluon, SCE for the scalar octet and do for the scalar singlet as 
indicated in Eq. ([2]). 





B{AA} 


B{AA,NS} 


■S{AA,iV3,SS} 


OGE 


R 


R 


R 


PSE 


R 


R 


R 




A(5.4) 


A(5.4) 


A(8.3) 


SCE 


WA 


WA 


R 


OGE + PSE 


R 


R. 


R 


do + PSE 


A(0.1) 


A(0.1) 


A(0.4) 


OGE + cj 


A(0.1) 


A(0.1) 


A(0.6) 


OGE + PSE + ao 


WA 


WA 


WA 


OGE + PSE + do + SCE 


WA 


A(1.6) 


A(7.0) 



different cases were considered: a pure one-gluon exchange model supplemented by the long 
range Yukawa part of the pseudoscalar interaction, a hybrid model containing both OGE 
and chiral pseudoscalar plus scalar interactions, and a pure chiral model. Among them, the 
pure one-gluon exchange model produced by far the largest binding energy. 



In our case the situation is qualitatively similar to Refs. [22J. The OGE interaction in the 
three diagonal channels is plotted in Fig. |2j The repulsive character of the gluon exchange in 
AA exceeds by far the attraction present in EE, and even more in NS. Therefore if one looks 
at the binding energy with only gluon interaction, one finds that the Fredholm determinant 
is larger than 1, what is associated to a repulsive force. This is a general conclusion valid 
for any value of a s , since it is a multiplicative constant that appears in the three channels 
and the ratio between the repulsion in AA and the attraction in EE is preserved. Indeed, 
such an approximated ratio {\Vqq E (R = 0)\/\Vqq E (R = 0)| ~ 5) is also found in (22|, where 
a SU(6) quark model is used, in spite of the huge difference in the quark-gluon coupling 
constant, taken there to be a s = 1.52, almost one order of magnitude larger. The global 
effect of the chromomagnetic interaction is thus found to be repulsive and therefore one has 
to keep looking for the responsible of the binding. 

In order to appreciate the overall features of any interaction, the S = —2, (T, S) = (0, 0) 
two-baryon system has been solved in the one- two- and three-channel approaches by 
including only the interactions indicated in each line of Table IIVI This gives us information 
on the character of the interaction involved, whether it is repulsive, weakly attractive or 
attractive enough to bind the system. The pseudoscalar interaction generates repulsion, 
and it contributes only to the AA channel, being negligible in the others. Therefore in a 
model containing only PSE, or OGE + PSE interactions there is not a bound state. The 



same situation occurs in Refs. 251 , where a resonance 26 MeV above threshold is found in a 



OGE + PSE model. However, the addition of a scalar exchange at the baryonic level provides 



large attraction so that the binding becomes Bu = 20 MeV. In Ref. [23[ the pseudoscalars 



also introduce repulsion at short distances, but weaker than in the NN interaction, in such a 
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0.5 0.52 0.54 0.56 0.58 0.6 



0.52 0.54 0.56 0.58 0.6 

g ch 2 (K)/4jt 



FIG. 3: Dependence of Bh on the quark-meson exchange couplings. See text for details. 



way that a bound H particle -provided there were medium range attraction as the a implies- 
is not discarded. Repulsive pseudoscalar contribution is also found in |24|, together with an 
attractive a piece, for a (Os) 6 configuration. When considering the two-cluster configuration 
with a pure chiral content the bound state does not appear. The reason for that is not a 
repulsive character of the chiral meson exchange but the lack of strength. Because of that, 
an increment of the chiral coupling constant in order to properly describe the 1 5'o NN phase 
shifts is successful for binding the system. 

In our model, the <7o is very attractive and provides by itself the largest values for Bh, 
from 5.4 MeV to 8.3 MeV. Let us notice that in a pure do model the two transition potentials 
to NE would not be possible, since they are due to the scalar meson exchanges. Then the 
AA would not be coupled to NE and so the two-channel calculation would give the same 
result as the pure AA, as can be checked in Table IIVI When adding the <7o to the OGE 
+ PSE terms, an unbound but near threshold H dibaryon is obtained in all cases. As in 
Ref. quite close to our treatment, the a Q is the only attractive piece. The contribution 
coming from the exchange of the scalar octet is weakly attractive or repulsive, as can be 
seen in Table IIVI However, when supplementing the OGE + PSE + ctq potential with the 
octet of scalars, a bound state appears in the two-channel calculation with B — 1.6 MeV 
and B = 7.0 MeV in the three-channel case. It is exactly at this point where the important 
role of the scalars in our model lies, or more precisely, that of the k exchange. This exchange 
piece manages to enlarge the binding energy without being attractive itself because it is the 
main piece in the AA-N3 and A^S-EE transitions. It redistributes part of the flux to the 
much more attractive NS and EE channels, originating thus an appreciable increase in the 
binding. Similar results have been found in the quark model framework of Ref. (26| . 



B. Parameter dependence 

As deduced from our discussion above, the binding energy of the dibaryon depends mostly 
on the o"o and k exchange potentials. Therefore it is interesting to evaluate the dependence 
of the binding on the parameters involved in these two interactions, more precisely on their 
masses and coupling constants. Varying the ctq and particularly the k coupling constants 
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FIG. 4: Dependence of Bh on the masses of (a) the a$ an d (b) k mesons. 



along a reasonable range changes the binding energy as seen in Fig. [31 being these changes 
larger in the first case. In both cases Bh grows with the coupling constant. Indeed it can 
be checked that the change on g c h, that affects the exchange of every meson, produces quite 
the same results as the change on g(o"o) exclusively. 

The dependence on the masses is opposite, see Fig. 0] for details. The binding decreases 
as the masses get larger, due to the shorter range of the attraction. This behavior is more 
pronounced in the case of m ao . Bh ranges approximately from 4 MeV to 12 MeV for 
m ao G (3.2, 3.6)fm -1 . Both dependence features (masses and couplings) are general, and 
similar qualitative results were also found in (24] |. 



C. Flavor symmetry breaking 

The flavor symmetry breaking (FSB) is expected to be an important effect in the H dibaryon 
due to the presence of two strange quarks. The present study has been made under the 
assumption that SU(3) is broken. Several sources of symmetry breaking are taking part in 
our work. First of all, the strange quark has been assigned a mass different to the light 
quark masses (m s = 555 MeV ^ m u ^ = 313 MeV). It enters in the expressions of the OGE 
and pseudoscalar potentials through the factor l/m^mj. However, the differentiation of the 
s-quark in the wave function through the harmonic oscillator parameter b s gives the larger 
source of FSB. Besides giving different values for the orbital matrix elements, it implies the 
suppression of a number of diagrams because light and strange quarks, as distinguishable 
particles, cannot be exchanged. 

In order to roughly evaluate the effect of FSB on the H dibaryon, the binding energy 
has been recomputed when some sources of symmetry breaking have been eliminated, i.e., 
for m s = m U ( i = 313 MeV and b s = b. We have variated in Fig. E]the harmonic oscillator 
parameter as well as the mass of the strange quark, between their values for exact SU(3) 
and their true values in our model. One can see how the binding energy increases as one 
approaches the exact SU(3) limit. As a consequence, the binding is increased by around 
45%, giving Bh = 10.2 MeV. Thus, we can conclude that FSB lowers the attraction. A 



similar conclusion was obtained in Ref. [20|, through a RGM calculation with only OGE 
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FIG. 5: Dependence of Bh on the strange quark (a) harmonic oscillator parameter and (b) mass. 



and confinement potentials. There, the binding energy drastically changed from 38 MeV for 
exact SU(3) to 26 MeV above threshold when the FSB effect was incorporated. A similar 



qualitative result was found by Ref. |24j, where the authors concluded that the symmetry 
breaking diminished the strong attraction coming from the chromomagnetic interaction. 
From chiral effective field theory j27|, also the bound state found for exact SU(3) disappears 
when using physical values for the baryon and meson masses. 

A closer inspection of every component of our potential in all the six channels contributing 
to H allows us to conclude that the effect of FSB is to diminish the strength of the potential, 
be it repulsive or attractive. Therefore, depending on the combination of the contributions 
in any channel, the total potential will be larger or smaller than the one with FSB, but 
no other general conclusion can be inferred. In our particular case, it turns out that AA 
and EE hardly change, whereas AS gains attraction when approaching SU(3). As for the 
transitions, the only appreciable change takes place in the EE-VH interaction, that acquires 
twice as much attraction for the exact SU(3) case. This result is in contrast with the one 
found in j22|, where the AA-VH interaction is strongly dependent on a FSB parameter. 



D. Spatial configuration 

At first sight there are two types of possible configurations for the H dibaryon: the six- 
quark cluster and the two-baryon cluster. As we have already mentioned the first one was 
employed in the original work [3|, getting a bound H. However latter calculations found that a 
dibaryon with only the (0s) 6 configuration would not bind [2l[. The spatial distribution was 
also investigated in Ref. [2J], where an extended resonating group method was employed 
in order to account for the possibility of a change in the baryon wave functions. Again, 
the dibaryon did not bind in the (Os) 6 configuration, being thus necessary to consider less 
compact configurations, due to the medium range attraction. When enlarging the size 
of the wave function, the dibaryon became more attractive, being Bh = 18.2 MeV for the 
stable solution. The probabilities however did not drastically change from the more compact 
((Os) 6 ) to the more spread (two-cluster) configuration. In the latter case they were found 
to be P A a = 0.216, P NS = 0.543 and P ss = 0.242. This result is slightly different from 
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FIG. 6: Interacting potentials contributing to the T = channels, (a) S = and (b) S = 1. 



that in the flavor singlet state, Eq. ffl~2]) . The small difference justifies the perturbative 
treatment employed, since the wave function was first postulated to be the flavor singlet 
when performing the RGM calculation. 

Usually, the flavor singlet wave function, Eq. ffl~2]) . is first postulated and a perturbative 
variational calculation is afterwards performed when making use of the RGM treatment. 
This is not the case of our work, since the coefficients of the baryon-baryon components 
of the flavor wave function are obtained as an output of the calculation, without making 
further assumptions. The probabilities we got are Paa = 0.177, Pnb = 0.446 and Pes = 
0.377. They are quantitatively similar to those of the flavor singlet, from what we can 
infer that in our model the baryon-baryon wave function is at first approximation SU(3) 
symmetric, being the difference between both due to the flavor symmetry breaking effects. 



E. Other states in S = — 2 

In order to properly understand the S = —2 sector, the analysis performed on the H dibaryon 
and the results obtained should be complemented studying the other spin-isospin channels. 
From such a complete analysis, more general and powerful conclusions may arise. 

Unlike the 1 Sq(T = 0) H dibaryon channel, the 3 S\(T = 0) is found to be very weakly 
attractive. Only the iVE interaction is present. The AE system does not couple to T = and 
the wave function of AA and EE is antisymmetric, and therefore vanishes, for L even. Since 
AA is forbidden, there are no lighter states and therefore the iVS state could decay weakly. 
The 3 Si(T = 0)NE potential, Fig. MJo), exhibits a repulsive core, due to the strong OGE 
interaction that compensates by far the attraction coming from the ao- At intermedium 
range the attractive feature survives but is not enough to form a bound state. The weak 
attraction in NS can be enhanced if the tensor force is turned on, however a bound state 



is still not possible. This qualitative behavior coincides with that of Refs. |22l . |28| and 
Nijmegen models D and F [29]. The big difference from the H channel is an example of the 
importance of the flavor dependence interaction that Model II incorporates. A summary of 
the predictions for the scattering lengths of the S = — 2 two-baryon systems with modern 
versions of the Nijmegen potential (30] . the Julich potential (HI], an d the work by Fujiwara 



et al. I26J can be found in Table 2 of Ref. 17 



The potentials in the two T = 1 channels are all attractive. Regarding the 1 So, both 
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FIG. 7: Interacting potentials contributing to the T = 1 channels, (a) S = and (b) S = 1. 



NE and AE interactions are strongly attractive, as can be seen in Fig. [7(a). Such a strong 
attraction provides for the coupled system a binding energy of 4.8 MeV. However this is 

23], where all the potentials with T = 1 are repulsive. As for the 

NE, AS and EE. The three of them are 
moderately attractive and look very similar, as can be seen in Fig. [7Jb). The Nijmegen 
models D and F 
in our model. 



not the case in Ref. 
3 5i channel, there are three coupled interactions 



301 ] also predict attractive interactions. However, no bound state appears 



IV. SUMMARY 



We have performed the first calculation of the H dibaryon in a model constrained by the 
elastic and inelastic AA^, UN, EN and AA cross sections. Special interest has been devoted 
to analyze the contribution of the different pieces of the interaction and to the effect that 
the addition of channels to the AA system produces on the binding. We obtained a bound H 
dibaryon, with Bh = 7 MeV, compatible with a plausible extrapolation of recent lattice QCD 
results and with the Nagara event, the most stringent restriction known so far. The scalar 
octet exchange, although not giving an attractive contribution by itself, plays a key role as 
it is the main ingredient of the transition potentials that connect AA to the more attractive 
NE and EE states. The probabilities P\a = 0.177, Pjvh = 0.446 and P^s = 0.377 are quite 
similar to those of the flavor singlet so that our wave function is at first approximation 
SU(3) symmetric. Finally, a bound state has also been found in the S = —2, (T, S) = (1, 0) 
channel, with a binding energy of 4.8 MeV, thus smaller than Bh- The abundance of events 
foreseen in a near future and/or the improvement on the lattice calculations will help us 
to advance in our knowledge of the mechanisms that play important roles in the dynamics 
of the H dibaryon. For such a task, a detailed theoretical study as the one presented here 
could be relevant. 
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